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On Infectious Model for Dependent Defaults 

Jia-Wen Gu * Wai-Ki Ching t Tak-Kuen Siu * Harry Zheng § 



Abstract 

In this paper, we propose a two-sector Markovian infectious model, which is an 
extension of Greenwood's model. The central idea of this model is that the causality 
of defaults of two sectors is in both direction, which enrich dependence dynamics. 
The Bayesian Information Criterion is adopted to compare the proposed model with 
the two-sector model in credit literature using the real data. We find that the newly 
proposed model is statistically better than the model in past literature. We also 
introduce two measures: CRES and CRVaR to give risk evaluation of our model. 

Keywords: Contagion Model, Markov Chain, Two-sector Model, Risk Management, 
Causality. 

1 Introduction 

Modeling dependent default risk has been a key issue in credit risk modeling. There 
are two important approaches to model the dependent default risk. The structural firm 
model has its origin in Merton (1974) and Black and Scholes (1973), which models the 
relationship between the firm's asset value and the defaults. The reduced-form intensity- 
based model by Jarrow and Turnbull (1995) use Poisson jump processes to model the 
default event. 
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Copula has been a very popular tool in modeling the dependent risk. The idea of 
Copula is transforming the marginal variables to uniform variables by a simple transfor- 
mation. After this is done, a n-dimensional function is used to model the dependence of 
the uniform variables, which is so called the Copula function. The Copula helps us to 
deal with the multivariate distribution of the uniform variable, without consideration of 
the original marginal variables. There are many useful Copulas in finance. The Gaussian 
Copula, which is introduced by Li (2000), is widely used in risk modeling and financial 
assessment. 

In addition, conditional independence model is also a commonly used model in credit 
risk modeling. Conditional on the systematical common factor, the loss random variables 
are independent. To specify, the Bernoulli mixture model is followed by the Credit Metrics 
and KMV-mode\, while the Poisson mixture model is followed by the CreditRisk + model. 
In a recession, the default of one company is triggered by the underlying common risk 
factor and also by the related company's defaults. The contagion model is used to de- 
scribe how the credit event of one company affects the other companies. Davis and Lo 
(2001) introduce an infectious default model, where in a portfolio a bond may be infected 
by defaults of other bonds or default directly. Jarrow and Yu (2001) propose a reduced- 
form model to describe the defaultable bonds of different company, where the concept of 
counterparty risk is first introduced to the credit literature. 

Ching et al. (2008) introduce an infectious default model based on the idea of Green- 
wood's model considered in Daley and Gani (1999) . This model aims at modeling the 
impact of default of a bond on the likelihood of defaults of other bonds. The original 
version of Greenwood's model is a one-sector model. It is then extended to a two-sector 
model in Ching et al. (2008). Besides, the joint probability distribution function for the 
duration of a default crisis, (ie, the default cycle), and the severity of defaults during the 
crisis period was also derived. Two concepts, namely, the Crisis Value-at-Risk (CRVaR) 
and the Crisis Expected Shortfall (CRES), are also used to assess the impact of a default 
crisis. The Greenwood's model is also extended to a network of sectors in Ching et al. 
(2010). Gu et al. (2011) propose a Markovian infectious model to describe the depen- 
dent relationship of default processes of credit securities based on Ching et al. (2008, 
2010), where the central idea is the concept of common shocks which is one of the major 
approaches to describe insurance risk. 

In this paper, we propose a two-sector Markovian infectious model, where the future 
default probability switching over time depends on the current number of defaults of both 
sectors. Moreover, the defaults of sector A caused by th defaults of sector B, and vice 
versa. The causality of defaults in both direction is captured by the underlying switched 
default probability. We adopt the maximum likelihood method to estimate the parameters 
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and the Bayesian Information Criterion to compare the propose model with two-sector 
model considered in Ching et al. (2008). The experiment result shows that the proposed 
model outperforms the model in credit literature. In addition, a more general model is 
given to provide more flexibility in describing realistic features of the dynamics of default 
probabilities. 

This paper is structured as follows. Section 2 presents our proposed model. And we 
also derive a recursive formula for the joint probability distribution for the default cycle 
and the number of defaults during the crisis and outline the estimation procedure. Section 
3 presents the ideas of the CRVaR and the CRES. In Section 4, we present the results of 
empirical analysis using our proposed model. Section 5 gives the general model extending 
the proposed model in Section 2. The final section concludes the paper. 



2 The Basic Model 

Let T be the time index set {0, 1,2,...,} of our model. To model the uncertainty, we 
consider a probability space (Q, J 7 , V), where V is a real-world probability. Suppose that 

X:={X t } t€T and Y := {Y t } teT 

denote two stochastic processes on (SI, J 7 , V), where X t = (X£,Xf) and Y t = (Y^,Yf) 
represent the numbers of surviving bonds and the defaulted bonds at t £ T in sector A 
and sector B, respectively, e.g., X\ represents the the number of surviving bonds at time 
t in sector A. We assume that the initial conditions are given as follows: 

X = (arj, Xq), Y = (y\, y 2 Q ) and x^ + yl = N u x 2 + y% = N 2 . 

Note that for each t £ 7", the sum of the numbers of the defaulted bonds and the surviving 
bonds at the time epoch t + 1 must equal the number of surviving bonds at time t in 
every sector, i.e., 

X] +1 + = X\ and Xl 1 + Y t \ 1 = Xl (1) 

For each t £ T, let at and flt be the probability that the default of a surviving bond is 
infected by the defaulted bonds at time t in sector A and sector B, respectively. The joint 
probability distribution of {X t+ i,Y t+1 } given {X t ,Y t } is given by the following Binomial 
probability: 

P(xt,y t )( x t+i,Vt+i) = P{(X t+1 ,Y t+1 ) = (x t+u y t+ i) | (X tl Y t ) = (x t ,y t )} 

= ( x } )k)^(i-^x( f (2) 

V vl+i j v y*+i J 
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We consider here the situation that the joint future default probability depends on the 
current number of defaulted bonds of both industrial sectors. We assume that 



a t 



a(Vt) 
a 
ai 
a 2 
a 3 



if 
if 
if 
if 



y] = yf = Q 
yl>0,y? = 
yt = 0,y?>0 
yl>0,y 2 t >0 



(3) 



aoh (yj, y\) + a x h x {y\, y 2 ) + a 2 h 2 (yj, y 2 ) + a 3 h 3 (yl, y} 



and 



Kvt) 

h 



if 
if 
if 
if 



y\ = yf = v 

Vl = 0,y?>0 
y}>0,y 2 = 

y\> 0,1/? >o 



(4) 



boh (y 2 t , y\) + hh^yl y\) + b 2 h 2 (y 2 t , y\) + 63^3(2/?, 2/t), 



where 



and 



ho{x,y) 



h 2 {x,y) 



if x = y = 
otherwise 

if x = 0, y > 
otherwise 



hi(x,y) 



h${x,y) 



if x > 0,y — 
otherwise 

if x > 0, y > 
otherwise. 



As it is shown in Equation ([I])©, one can see that {X t ,t = 0, 1, 2, . . .} is a second- 
order Markov chain process. We remark that this two-sector model provides a novel 
and flexible dependent structure for correlated defaults of two different industrial sectors. 
Firstly, an infectious default within one time period is modeled as a Binomial distribution, 
which has been widely used in modeling the spread of epidemics whose situation is quite 
similar to that of a financial crisis. The causality of the infection is supposed to be in 
both direction, i.e., a "looping default". Secondly, the process (X t ,Y t ) has the Markov 
property, where the probabistic structure of future states only depend on the current 
state. Thirdly, conditioning on the current state (X t ,Y t ), the future state of two sectors 
(Xh^Y^) and (X 2 +1 ,Y 2 +1 ) are stochastically independent. The step functions hi(x,y) 
are used to describe the dependence of the default probabilities on the state of previous 
time epoch. On one hand, this method provides a tractable and analytic solution for 
parameter estimation from empirical data. On the other hand, one has to admit that this 
simplicity may result in limitations in applications. In Section 5, we relax the assumption 
of the specific form for at and fit and a more complicated dependent structure modelling 
framework is presented. 
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2.1 Default Cycle and Severity 

In this subsection, we proceed to derive the joint probability distribution function (p.d.f ) 
for the duration of the default crisis (T), namely, the default cycle, and the severity of the 
defaults (Wt) during the crisis period. These two concepts are essential in determining 
the impact of a default crisis. We first give a precise definition of the default cycle: 

T := mf{t e T | Y t = 0}. (5) 

And given T = t > 0, W t represents the number of defaults in the sector over the time 
duration (0, t]. To apply the concepts of default cycle and the severity of the defaults on 
our proposed two-sector model, we write 

T\ := inf{* e T | Y t l = 0} and T 2 : = inf{t eT\Y t 2 = 0}. 

Provided that T\ — t± > and T 2 = t 2 > 0, Wf and W? represent the number of defaults 
in sector A and sector B respectively in (0, ti] and (0, £2]- To obtain the joint distribution 
of (Wj. , Ti) for % — 1, 2, we assume that (X , Y ) = (x , y ) with > 0, y\ > 0. Let 

P n (xi,x 2 ,h) = P{T X >n + l,X* = xi,X* = x 2 ,I {Y 2>o} = h}. 

The following Lemma gives recursive formulas for P n (xi,X2, h). 

Lemma 1 



P n (x u x 2 ,0)= J2 ( Sl ) \Pu^(suX 2 ,0)(a i y^(l-a 1 )^(l-b 2 ) 



Sl>Xl 



■r-2 



J-2 



+P n -i( Sl , x 2 , l)(a 3 ) s ^(l - a 3 r(l - b 3 ) 
P n (x u x 2 ,l)= J2 E f Sl V S2 V P -i( s i' s 2,0)(a 1 ) Sl -^(l-a 1 )^(& 2 ) S2 - a;2 (l-&2 

81>X1 Sl >xi \ X lJ \ X 2/ 

+p n -x{s 1 , 8 2 , i)(a 3 y^(i - a 3 )^(b 3 y^(i - b 3 y> 



I .('2 



where the initial condition is given by 

P (x 1 ,x 2 ,h) = 



1, (x u x 2 ,h) = (x}„xl, 1) 
0, otherwise. 
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Proof: By the law of total probability and Markov property, 
P n (x 1 ,x 2 ,0) 

= P{T X >n+l,X 1 n = x 1 ,Xl = x 2 , 1 {Y 2 >0} = 0} 

= E E P { T i ^ n > X n-i = si, X 2 n _ x = x 2 , I{Y*_ t >o} = h} 

s 1 >xi h=0,l 

xP{Ti > n + 1, Xl = xi, Xl = x 2 , I{y%>o} = | Ti > n, Xl_ x = s 1} X$_ x = x 2 , Zp^x)} = M 
= E E P n -i(si,x 2 ,h) 

si>xi h=0,l 

xP{Y^ > 0, X\ = xi,X% = x 2 , /{y n 2>o} = | 71 > n, X\_ x = Sx,Xl_ x = x 2 , I{ Y ^_ 1>0 } = h} 
= E E P n -i(si,x 2 ,h) 

s 1 >xi h=0,l 

P{Y^ > 0,Xl = Xl , XI = X 2 ,I {Y 2 >0} = | Y*-x > = S l>*n-1 = X2j{Y^_ 1>0 } = h} 



= E (: i 1 )[^i(5i,x 2 ,0)(a 1 ) s -^(l-a 1 )^(l-fo 2 

Sl>X'l 

+P n _x(sx,x 2 , l)(a 3 ) Sl ^(l - a 3 ) xi (l - h) x * 
Similarly, we have 
P n (xx,x 2 , 1) 

= P{T 1 >n + l,X 1 n = x 1 ,X? l = x 2 ,I {Y 2 >0} = l} 

= E E E P { T l>^^-l = S l^n-l = S 2,/ { y 2 _ 1 >0} = M 
si>a;i S2>Z2 h=0,l 

xP{Ti > n + 1, X^ = xi, Xl = x 2 , /{y2> } = 1 | Tx > n, X\_ x = s u Xl_ x = s 2 , /{y2_ 1>0 } = h} 

= E E E Pn-l(sx,S 2 ,h) 
si>xi S2>X2 h=0,l 

xP{Y' >0,X* = xx, Xl = x 2 ,I {Y 2 >0} = l\Tx> n,Xl_ x = Sl ,Xl_ x = s 2 ,I {Y *_ i>0} = h} 

= E E E Pn-l(sx,S 2 ,h) 
si>X! s 2 >x 2 h=0,l 

P{Yl > 0,Xl = Xl ,Xl = x 2 ,I {Y 2 >0} = 1 | Y*_i > 0,Xl_ x = sx,Xl^ = s 2 ,I {Y 2_ i>0} = h} 



E E (:;)(: 2 J[Pn-i( S i, S2 ,o)(a 1 )^Mi-air i (^r 2 --(i-^; 

Sl>Xl Sl>Xl 

+P n -x(sx, s 2 , l)(a 3 r-^(l - a 3 )^(b 3 ) s ^(l - b 3 f 



□ 



Proposition 1 The joint distribution of (Tx,Wj. ) is given by 

P{(Tx, W^) = (n,x)} = Y J Pn-i{xl-x, x 2 , 0)(l-ax) x °- x +Yl Pn-i(4~^ x 2 , l)(l-a 3 ) x °- 

X2 X 2 



6 



Proof: 

P{(T u Wl) = (n,x)} 
= P{T 1 >n^=Q,X l n = xl-x} 

= P { T i ^ n > X n-i =x\-x, X 2 n _ x = x 2 , ^{y n 2_ 1>0 } = h} 

X2 h=0,l 

x P{Yi = 0,X^ = x 1 o -x\T 1 >n, X l n _ x =x\-x, X$_ x = x 2 , I {Y ^_ 1>0} = h} 

= H Pn-l{xl -X,X 2 ,h) 

X2 h=0,l 

*P{Yn =0,X* = xl-x \ Y^_ x > 0, X l n _ x = x 1 -x, X\_ x = x 2l Ip*_ l>0} = h} 
= Ex 2 P n -\{xl - x, x 2 , 0)(1 - a x )<- x + J2 X2 Pn-i{xl - x, x 2 , 1)(1 - a 3 ) x °- x 

□ 

We remark that due to the symmetric property of the two sectors, the joint distribution 
(M, T 2 ) shares a similar form of (W^ , T\). 



2.2 Parameter Estimation 

This two-sector model has eight parameters: a , a±, a 2 , a 3 and b , b x , b 2 , b 3 . We 
employ the maximum likelihood method to estimate the parameters. Given the total 
bonds Ni, N 2 and the observations of the number of defaulted bonds y\, y\, . . . , y l N and 
y%, yf, . . . , yff, where N denotes the period of observation time, the number of surviving 
binds x , x\, . . . , x N and deterministic. The following proposition gives 

analytical expressions for the maximum likelihood estimates of the model parameters. 

Proposition 2 For i — 0, 1, 2, 3, 

JV-l N-l 

- - 4=0 and bi - 4=0 



JV-l 1 iV-l 

E xlWvly}) £ x 2 Mylyl) 

t=o t=o 

Proof: We prove the expression for So here and the proof for the others are similar. 

The likelihood function L(a, b \ x , x±, . . . , x^, y , y±, . . . , y^) is then the joint probability 

density function f(x , xt,...,x N , y , y u . . . ,y N \ a,b): 

L(a, b | x ,X!, . . . , x N ,y , y u . . . ,y N ) 
= f(x , x 1 ,...,x N ,y ,yi,...,y N \a,b) 

a(y ))MyoY l x [ i ) (1 - b(y )rU(y f 





a( yi ))Myif 2 x [ M (i - b( yi )rlb( yi f* 

- a( 2/jv _ 1 ))^a(y JV _ 1 )^ x ( aJ *" 1 ] (1 - b( 2/Ar _ 1 ))^b( 2 / JV _ 1 )^. 



N 
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Then by solving 

din L(a,b | x ,x l7 . . .,x N ,y ,y u ...,y N ) _ 
<9an 

we have 

& S$ a(2/t) 

Since for any t, 

1 _j^ h i (ylvl) and 1 _y h{y},yt) 
1 - a (?/t) i^n 1 ~ a * a (^*) ^ ai 



then 



Thus, 



i=o x ^\yt/ i=o 



N £j, xlMvlv?)hi{yhvl) , N ^j^ ylMyly 2 t)h,{yhy 2 t ) 

1 - an 



a t=n a 



t=0 x t=n 



AT-l 

E yl+Myhy 2 ) 
E xtho(yl,yi) 

t=0 



□ 



3 Crisis VaR and Crisis ES 

In this section, we give a brief introduction to the concepts of the CRVaR and the CRES 
in Ching et al. (2010). Then we present the evaluation of the CRVaR and the CRES 
using the proposed models. The CRVaR and the CRES are measures for the duration 
and the severity of a default crisis. Let 

L(-,-)(u) :T x K x n -> TZ 

be a real- valued function L(T,Wt)(uj) of T and Wt- We then suppose that for a fixed 
uj G f2, 

r(cd) = £, W t (u;) = w, and L(t, w)(u) = l(t, w) G TZ. 

That is, the loss from the default crisis is /(£, w) when the duration of default crisis T = t 
and the number of defaulted bonds in the crisis W t = w. We write jC(T, Wt) for the space 
of all loss functions L(T, Wt){oj) generated by T and Wt- 
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Sectors 


lotal 


Defaults 


Consumer 


1041 


251 


Energy 


420 


71 


Media 


650 


133 


Transport 


281 


59 



Table 1: The default data (Taken from Giampieri et al. (2005)). 

The CRVaR with probability level (3 under V is then defined as a functional V^(-) : 
C(T, W T ) ->> n such that for each L(T, W T ) G C(T, W T ), 

Vp(L(T, W T )) ■= inf{Z G 1Z\V(L(T, W T ) > I) < P} . (6) 

In the language of statistics, Vp(L(T, Wt)) is the generalized /3-quantile of the distribution 
of the loss variable L(T, Wt) under V. Since the loss from the default crisis L(T, Wt) 
is completely determined when T and Wt are given, V(L(T, Wt) > I) is completely 
determined by the joint p.d.f. of Wt and T. 

The CRES with probability level /3 under V is also defined as a functional -Eg(-) : 
C(T, W T ) -> U such that for each L(T, W T ) G C(T, W T ), 

Ep(L(T,W T )) := Ev[L(T,Wt)\L(T,W t ) > Vp(L(T,W T ))\. (7) 

In other words, Ep(L(T, Wt)) is the average of the loss from the default crisis when the 
loss exceeds the CRVaR of the default crisis with probability level f3 under V . 

4 Empirical Results for Proposed Model 

In this section we present the empirical results of the proposed two-sector model using 
real default data extracted from the figures in Giampieri et al. (2005), where we adopt 
the estimation methods and techniques presented in the previous section. 

The default data comes from four different sectors. They include consumer /service 
sector, energy and natural resources sector, leisure time/media sector and transportation 
sector. Table 1 shows the default data taken from Giampieri et al. (2005). From the 
table, the proportions of defaults for Consumer, Energy, Media and Transport are 24.1%, 
16.9%, 20.5% and 21.0%, respectively. The default probabilities of all four sectors are 
significantly greater than zero. This means that the default risk of each of the four 
sectors is substantial. 

We then construct the infectious disease model using these real data. The asterisk "*" 
in the table indicates the pair of sectors which has the largest correlation. From Table 
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Energy Consumer 




Transport Media 



Figure 1: The partner relations among the sectors using correlation. 





Consumer 


Energy 


Media 


Transport 


Consumer 




0.0224 


0.6013* 


0.3487 


Energy 


0.0224 




0.1258* 


0.1045 


Media 


0.6013* 


0.1258 




0.3708 


Transport 


0.3487 


0.1045 


0.3708* 





Table 2: Correlations of the sectors. 

2, we see that all correlations are positive. This provides some preliminary evidence for 
supporting the use of the two-sector model from the perspective of descriptive statistical 
analysis. We shall provide more empirical evidence for supporting the use of the proposed 
infectious model by the results of Bayesian Information Criterion (BIC) later in this 
section. To build the infectious model, for each row (sector A), we may find a partner 
(sector B) by searching the one with the largest correlation in magnitude (ie, the one with 
the asterisk "*"). Figure 1 gives the partner relations among the sectors using correlation. 
Later in this section, we will give the results for BIC to support the matched pair presented 
in figure 1. The estimation results for proposed infectious model and two-sector model 
Ching et al. (2010) are presented in Table 3. 

To compare the proposed infectious model with the two-sector model Ching et al. 
(2010), we consider the Bayesian information criterion (BIC). The formula for the BIC is 

BIC = -2log(L) + klog(n), 

where n is the number of observation data, k is the number of free parameters to be 
estimated, and L is the maximized value of the likelihood function for the estimated 
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Sector A 


Consumer 


Energy 


Media 


Transport 


Sector B 


Media 


Media 


Consumer 


Media 


Proposed Model 


a 


0.0007 


0.0004 


0.0005 


0.0013 


a x 


0.0018 


0.0033 


0.0005 


0.0012 




0.0013 


0.0018 


0.0017 


0.0026 




0.0049 


0.0032 


0.0042 


0.0052 


Two-sector Model [2] 


a 


0.0013 


0.0018 


0.0005 


0.0013 




0.0043 


0.0023 


0.0033 


0.0036 



Table 3: Estimation Results for Proposed Model 



model. Given any two estimated models, the model with the lower value of BIC is the 
one to be preferred. Table 4 presents the value of the BIC for the proposed model and 
the two-sector Ching et al. (2010). We remark that for all the four sectors, the proposed 
model with lower value of BIC is statistically better. 



Sector A 


Consumer 


Energy 


Media 


Transport 


Sector B 


Media 


Media 


Consumer 


Media 


BIC (proposed model) 


419.0813 


215.4654 


301.2534 


2.1287 


BIC (two-sector model Ching et al. (2010)) 


434.6700 


231.8225 


321.0501 


2.1460 



Table 4: The Value of BIC for Proposed Model and Two-sector Model Ching et al. (2010) 



To compare the matched pairs in Figure 1 with other matched pairs for the proposed 
model, we also adopt the BIC. Since the models of different matched pairs have the 
same number of parameters and length of data set, to compare their BIC is equivalent 
to compare their log-likelihood ratio. Table 4 presents the log-likelihood ratios for the 
matched pairs in Figure 1 against other matched pairs. We remark that all the log- 
likelihood ratios are positive which support the matched pairs in Figure 1 for the proposed 
model. 

Our proposed model aims at modeling causality of defaults in both direction. From 
the pair up results, one may found that the relation is not necessarily symmetric. This 
relation is only found symmetric for the sectors media and consumer, which means the 
causality of defaults from both direction is more reasonable for the media and consumer 
sector. 

We provide a scatter plot to depict the correlation of defaults in the matched sectors. 
A simulation of defaults in matched sectors in our proposed model is also conducted. 
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Matched Pairs in Figure 1 


Sector A 


Consumer 


Energy 


Media 


Transport 


Sector B 


Media 


Media 


Consumer 


Media 


Other Matched Pairs 


Sector A 


Consumer 


Energy 


Media 


Transport 


Sector B 


Energy 


Consumer 


Energy 


Consumer 


log-likelihood ratio 


33.1330 


7.3286 


18.6264 


1.9942 


Sector A 


Consumer 


Energy 


Media 


Transport 


Sector B 


Transport 


Transport 


Transport 


Energy 


log-likelihood ratio 


10.7231 


7.3495 


14.6136 


8.4934 



Table 5: The Value of BIC for Matched pairs in Figure 1 and Other Matched Pairs 



Figure 2 presents the number of surviving bonds in the matched sectors of empirical data 
and simulation. 

To apply the two measures CRVaR and CRES in the proposed model, we consider 
some hypothetical values for the loss. The loss L(Wt, T), for each T = 1, 2, . . . , X and 
Wt = 0, 1, . . . , Xq, are as in flS]). Then we present the value of CRVaR and CRES for the 
proposed model as well as the two-sector model Ching et al. (2010) in Table 6. And the 
loss distribution are presented in figure 3. 

L(0,j) = j -1 + 0.1, for each j = 1,2,..., X ; 
L(i,j) = L(0,j) +2 — 1, for each i = 1, 2, . . . ,Xq and j = 1, . . . , Xq. 

From Table 6, we see that for all of the four sectors, the existing two-sector model 
underestimates both the CRES and CRVaR. This reflects that failure to incorporate the 
contagion effect described in our proposed model leads to an underestimation of credit risk 
and has important consequences for credit risk management, such as inadequate capital 
charges for credit portfolios. Indeed, the loss distribution implied by the proposed model 
has a much fatter tail than that arising from the existing two-sector model. This explains 
why the proposed model provides more prudent estimates for the risk measures than the 
existing two sector model via incorporating contagion. We also remark that the contagion 
model including the causality of defaults in both direction, (i.e., looping defaults), has a 
significant impact on the loss distribution. 
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Figure 2: Number of Surviving Bonds in Matched Sectors. 

5 A Generalized Model 

As in the basic model, the stochastic process (X t ,Y t ) has the Markov property, where 
conditioning on (X t , Y t ), (X^ +1 , Y^ +l ) and (Xf +1 , Y^ +1 ) are stochastically independent. The 
joint probability distribution, given the realization of (X t , Y t ), a t , /3 t , is given by: 

P( xt ,y t ),( at ,f}t)(xt+i, Vt+i) = P{{X t+1 , Y t+1 ) = (x t+1 , y t+1 ) | (X t , Y t ) = (x t , y t ),a t , /3 t } 

= f i ]K)^(l-a 4 )^x f f ) (P t )Ml 

However, instead of maintaining the specific form of a bivariate step function for at,(3t, 
this model assumes that a t and (3 t follow certain Beta distributions depending on (X t , Y t ). 
By assuming a beta density on the unknown transition parameters, the chain becomes 
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Sector A 


Consumer 


Energy 


Media 


Transport 


Sector B 


Media 


Media 


Consumer 


Media 


Proposed Model 


CRVaR(/3 = 0.05) 


374.1 


25.1 


122.1 


26.1 


CRES(/3 = 0.05) 


424.7 


33.8 


150.4 


33.8 


CRVaR(/3 = 0.01) 


457.1 


39.1 


168.1 


39.1 


CRES(/3 = 0.01) 


495.1 


47.5 


192.4 


46.5 


Two-sector Model Ching et al. (2010) 


CRVaR(,5 = 0.05) 


114.1 


12.1 


34.1 


10.10 


CRES(,5 = 0.05) 


146.1 


17.1 


45.7 


14.1 


CRVaR(/3 = 0.01) 


166.1 


20.1 


52.1 


16.1 


CRES(/3 = 0.01) 


195.6 


24.5 


63.3 


20.2 



Table 6: CRVaR and CRES 



a Markov chain with transition matrix containing random parameters. This allow us 
to incorporate parameter uncertainty while, at the same time, retaining the analytical 
tractability of the model. For each time period, the number of defaults has the Beta- 
binomial distribution depending on the number of defaults in last time period. The Beta- 
binomial distribution is extensively used in Bayesian statistics, empirical Bayes methods 
and classical statistics as an overdispersed binomial distribution. 

Specifically, it is assumed that the density of a t and j3 t are given by f a (x; (X t , Yt)) and 
fp(x] (X t ,Y t )), respectively, and 



Hvhvi 



2-) 



f a (x; (X t ,Y t )) = £ 



X 



and 



where 



i=0 B \ B H^ B i2) 



B(x,y) = ( f-Hl-ty^db 
Jo 



and Aij,Bij,i = 0, 1,2,3, j = 1,2 are parameters of the Beta distribution. From the 
definition, one can have the following transition probability: 
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;ure 3: Loss distribution for proposed model and two-sector model Ching et al. (2010). 



P{{Xl +1 X + i) = 04h,vW I (Xt,Y t ) = (x t ,y t )} 
' E 



vl+i 



. 1 L 



(a t yhi(l- at yt + i\(X t ,Y t ) = (x t ,y t ) 
p yl t+i(l-p) xl *+if a (p; (x t ,y t ))dp 



(10) 



vl+i ) t^ B i A iu A i2) Jo 



\A i2 -l 



dp 



hu,vi) B(v H^ x k +Am) 



Vt+i j i=o B(A il ,A l2) 
A similar transition probability distribution is shared with the number of defaults in 
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Sector 13 

PiiXl^li) = (4fi^+i) I (*t,Y t ) = (x t ,y t )} 

A U, /2 1 ^(^ 2 + i + ^i,^ + i + ^ 2 ) (ii) 

Z/t+i / i=o B(B a ,B i2 ) 
From the transition probability distribution, to obtain a closed-form solution for a max- 
imum likelihood estimate is difficult, if not impossible. However, we can compute the 
maximum likelihood estimates using numerical optimization. 

5.1 Default Cycle and Severity 

To derive the joint distribution of (WJ.,Tj) for i = 1,2, we repeat the same steps in 
Section 2 to compute 

P n (x 1 ,x 2 ,h) = P{71 >n+l,X* = Xi,Xl = x 2 ,I {Y 2 >0} = M- 



Lemma 2 

P n (x 1 ,x 2 ,0) 



Sl 



si>a;i \ X 1 



P n _i(si,a;2,0) 
P n -i(s 1 ,x 2 , 1 



B(si -xi + An, xi + A 12 ) B(B 21 ,x 2 + B 22 ) 



B(A n ,A 12 ) B{B 21 ,B 22 ) 
B( Sl - x x + A 31 , xi + A 32 ) B(B 31 ,x 2 + B 32 ) 



B(A 31 ,A 32 ) 



B(B 3 i, B 32 ) 



(12) 



PjXl,X 2 , 1; 



^ In I To) n ' l[ U 2 ' j B(Aii,Aia) B(B 21 ,B 22 ) 

»i>xi si>xi \ x 2/ 



-P n _i(si, S 2 , 1) 



B(s 1 -x x + A 31 , xi + A 32 ) B(s 2 -x 2 + B 31 ,x 2 + B 32 ) 



B(A 3U A 



32 J 



B(B 31 , B 32 ] 



(13) 



where the initial condition is given by 

P (x 1 ,x 2 ,h) = 



1, (xi,x 2 ,h) = (4,^0,1) 
0, otherwise 



Proof: We prove the first equality. The proof of the second one is similar. As in the 
proof of Lemma 1, 

P n (x 1 ,x 2 ,0) 

= J2 Pn-\{s\iX 2l h) 

si>xi h=0,l 

xP{Y^ > 0,X* = x x , Xl = x 2 ,I {Y 2 >0} = | Yl_x > 0,X^ = s^X^ = x 2 ,I {Y 2_ i>0} 
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Note that by (ITU1) and ([Til, 

P{Y* >0,Xl = x u X^ = x 2 , 1 {Y 2 >0} = | Yl_ x > 0, Xl_ x = s u Xl_ x = x 2 , / { ^_ 1>0 } 
fsA B( Sl -xi + A u , xi + A 12 ) B(B 21 ,x 2 + B 22 ) 
W B(A n ,A 12 ) B(B 21 ,B 22 ) 

and 

P{Yn > 0, X* = x 1; X^ = x 2 , 1 {Y 2 >0} = | > 0, X\_ x = snX^ = x 2 , 1 {Y ^ 1>0 } 
f Sl \ B( Sl - zi + A 3 i, xi + A 32 ) B(B 3U x 2 + ff 32 ) 
W S(A 3 i,A 32 ) B(B 31 ,B 32 ) 

Combining these two results, ( fl~2l) follows. □ 
Hence the joint distribution of (Ti, VF^J follows 

Proposition 3 

B(A u ,xl -x + A 12 ) 



P{(T 1; M^) = (n, x))} = J2 P n-i(4 ~x,x 2 , 0) 



x 2 B(A n ,A 12 
+ ^P n _iOo - x,x 2 , 1) 



S(A 31 ,4-x + A 32 ^ 



(14) 



B(A 31 ,A 32 ) 
Proof: As in the proof of Proposition 1, 



P{(T u Wl) = (n,x)} 

= S P n-l(xl ~ X,X 2 ,h) 

X2 h=0,l 

xP { Y n = 0, Xl = xl - x | Y^_i > 0, Xl_ Y = x\-x, Xl_ x = x 2 , / { y n 2_ 1>0} = &} 
Note that by (1T0I) , 

P^ 1 = 0, X* = ^ - x | y^x > 0, X^ = x* - x, X^ = x 2 , / { y 2 _ i>0 } = 0} 
B(A 11 ,^-x + A 12 ) 
5(^ U ,A 12 ) 

and 

P{Y^ = 0,X' = xl-x\ Y^ > 0,X^ = xj - x,X^ = x 2 ,/ { y 2 _ i>0} = 1} 
P(A 31 ,xi - x + A 32 ) 

5(^31,^2) 

Combining these two, ( fT4l) follows. □ 



17 



6 Concluding Remarks 



We propose a two-sector Markovian infectious model. The proposed model incorporated 
two important features of credit contagion, namely, the chain reactions of defaults and 
the bi-lateral causality of defaults between two industrial sectors. We capture the chain 
reactions of defaults by postulating that the future default probability switches over time 
according to the current number of defaults of two industrial sectors. The bi-lateral causal- 
ity of defaults meant that defaults in one sector are caused by defaults in another sector, 
and vice versa. This bi-lateral causality of defaults enriches the dependent structures 
of credit risk model. We provide an efficient estimation method of the proposed model 
based on the maximum likelihood estimation. Two important risk measures, namely, the 
CRVaR and the CRES, were evaluated under the proposed model. To provide a more 
flexible and realistic modeling framework for the dynamics of default probabilities, we 
extend the model to a case where default probabilities are Beta random variables given 
the realization of the state in the previous time period. 

We also conduct empirical studies on the credit risk models using real default data. 
We adopted the BIC to compare the proposed model with the existing two-sector model 
proposed in Ching et al. (2010). The numerical results reveal that the proposed two- 
sector model outperforms empirically the existing model. By comparing the risk measures 
evaluated from the proposed model and those evaluated from the existing two-sector 
model, we found that failure to incorporate the contagion effect described in the proposed 
model leads to an underestimation of risk measures. This provides some evidence to 
support the proposed model. 

One possible topic for future research may be to incorporate the impact of the number 
of defaults on the likelihood of future defaults via a different parametrization of the future 
default probability. In current paper, we assumed that the joint future default probability 
switches over time depending on the region where the current number of defaults falls in. 
Four parameters, namely, a , a 1; a 2 and were involved. To provide a more parsimonious 
way to incorporate the current number of defaults on the joint future default probability, 
one may consider the following parametrization for the future default probability: 

a t = a + aiy l t + a 2Vt > 

where y\ and y\ are the current numbers of defaults in the two industrial sectors. Using 
this parametrization, we can reduce the number of parameters by one and accounts for 
more information of the current number of defaults when evaluating the future default 
probability. 
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